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The chiral phase transition at finite temperature is investigated in the linear sigma model, 
which is regarded as a low energy effective theory of QCD with three momentum cutoff, 
in the variational method with the Gaussian approximation in the functional Schrodinger 
picture. It is shown that the Goldstone theorem is retained and the meson pair excitations 
are automatically included by taking into account the linear response to the external fields. 
It is pointed out that the behavior of chiral phase transition depends on the three-momentum 
cutoff, which leads to the careful treatment of the problem. 
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Q_i' §1. Introduction 

D \ It is well known that the chiral symmetry plays an essential role in the quark 

, ' and hadronic world governed by the quantum chromodynamics (QCD). Especially, 

the spontaneous chiral symmetry breaking leads to the existence of the pion, 1 ) which 
plays an important role in the nuclear and hadron physics. One of the recent interests 
in the hadron physics is the chiral symmetry restoration in the context of the BNL 
Relativistic Heavy Ion Collider (RHIC) and/or the CERN Large Hadron Collider 
(LHC). 2 ) It is believed that the chiral symmetry is restored in the high temperature 
and/or high baryon density, which might or may be created in the above mentioned 
experiments. 

00 \ The finite temperature phase transition in the field theoretical treatment was 

investigated in the context of the phase transition in the early universe. 3 ) System- 
atically, the field theoretical treatment for the systems at finite temperature was 
developed in the early stage. 4 )' 5 ) Further, another approach 6 ) to the finite temper- 
ature systems was devised with the variational technique in terms of the functional 
Schrodinger picture with Gaussian approximation 7 )' 8 ) for the bosonic field theory. In 
this paper, we intend to adopt a slight different variational approach to the bosonic 
field theory at finite temperature. 

As for the low energy region of the strong interaction, the low energy effective 
theories of QCD give very useful tool until now to understand the physics of hadrons, 
especially, based on the chiral symmetry, while the lattice QCD simulation is devel- 
oped for the problems related to the chiral symmetry. In the low energy region of 
QCD, the pion and its chiral partner, namely sigma meson, are important ingredi- 
ents. As an effective theory of QCD, the linear sigma model 9 ) gives a powerful tool 
to understand the many properties of low energy hadron physics, which contains 
pions and sigma meson paying an attention to the chiral symmetry. 

In the linear sigma model, the chiral phase transition at finite temperature was 
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firstly investigated in Ref.10), in which approximations developed in the many-body 
theory were used widely. However, according to the approximations, the order of 
the chiral phase transition was not settled definitely. Recently, in the context of 
the relativistic high energy heavy ion collision experiments, the renewed interests 
are pointed to the chiral phase transition in the O(N) linear sigma model from 
the theoretical sides. 11 )' 12 ) In the variational approach to the O(N) linear sigma 
model, it has been pointed out that the Goldstone theorem 13 ) is not retained in 
naive treatments. For solving this problem, some works have been carried out. 14 ) -17 -* 
Also, one of the present authors (Y.T.) with Vautherin and Matsui gives another 
method to recover the Goldstone theorem in the variational approach to the O(N) 
linear sigma model, in which to take into account the linear response for the external 
fields is essential. 18 ) 

As for the subject of the chiral phase transition at finite temperature, in Ref.ll), 
so-called optimized perturbation theory is constructed for the O(N) linear sigma 
model, in which a certain resummation of the higher order terms is carried out. In 
the case N = 4 for the realistic world, the Goldstone theorem is satisfied in this 
approach and it has been shown that this approach has revealed the second order 
chiral phase transition at finite temperature with realistic pion mass. On the other 
hand, the first order phase transition occurs with low pion mass including the case 
of the chiral limit. Another approach to the same subject was given in Ref.12) which 
gave a self-consistent Hartree approximation for meson mass equations for one-loop 
effective potential. Both approaches are carried out the renormalization procedure 
because the O(N) linear sigma model is renormalizable. 

As another approach, the large N limit gives interesting features for O(N) linear 
sigma model. At zero temperature, the vacuum stability was first investigated. 19 ) For 
large N limit in the O(N) linear sigma model, the Goldstone theorem is not broken. 
In the context of the chiral restoration, the large N approach was investigated and the 
order of phase transition is discussed 20 ) by using of the Cornwall-Jackiw-Tomboulis 
(CJT) effective potential. 21 ) 

In this paper, we investigate the chiral phase transition at finite temperature 
based on the time-dependent variational approach with Gaussian wave functional 
and with the imaginary time formalism of Green's functions, 22 ) in which the linear 
responses for the external fields are taken into account. At zero temperature, it has 
already been shown that the Goldstone theorem is retained in our approach for the 
O(N) linear sigma model with N = 4 mentioned above. 18 ) In this paper, we extend 
this approach to the finite temperature case and show the behavior of the chiral 
phase transition. We regard the 0(4) linear sigma model as a low energy effective 
model of QCD, so we introduce the three momentum cutoff parameter A which cuts 
the ultraviolet divergences of the meson loop appearing in our approach. It is shown 
that the pion-pion, pion-sigma meson and sigma-sigma mesons pair excitations are 
automatically included in our approach, which results to the Goldstone theorem and 
leads to the second order chiral phase transition with a reasonable parameterization 
in the 0(4) linear sigma model. 

This paper is organized as follows: In the next section, we recapitulate the basic 
ingredients of the variational approach with the linear response for the external fields 
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to the O(N) linear sigma model. The pion and the sigma meson masses are given 
and it is shown that the meson pair excitations are automatically taken into account. 
In §3, we extend this approach to the finite temperature system in the Matsubara 
formalism. 22 ^ In §4, some numerical results and discussion are given. The last section 
is devoted to a summary and concluding remarks. 

§2. Variational approach to the 0(4) linear sigma model satisfying the 

Goldstone theorem 

In this section, we recapitulate the basic ingredients of the variational method 
to the O(N) linear sigma model with the Gaussian wave functional and with the 
polarization tensors appearing in the linear response theory for the external fields. 

2.1. Variational approach and the breaking of the Goldstone theorem 

We summarize the time-dependent variational method following the Refs.24) 
and 18). 

Let us start with the following Hamiltonian of the 0(4) linear sigma model with 
the explicit symmetry breaking term: 

H(x) = H (x) - cip a {x)5 a0 , 

11 1 A 

H (x) = -■n a (x)-K a (x) + -Vip a (x)-Vip a (x) + -rr^(pa(x)(p a (x) + — (ip a (x)(p a (x)) 2 . 

(2-1) 

Here, a runs from to 3, where ipo means the sigma field operator and ipi with 
i = 1 ~ 3 mean the pion field operators. In this paper, we adopt the time-dependent 
variational approach in the functional Schrodinger picture. In this picture, the conju- 
gate operators ir a (x) are replaced by the functional derivative ir a (x) = —i6/6ip a (x). 
A possible trial wave functional including the quantum fluctuations is a Gaussian 
wave functional which is written as 

¥[tp] = (<p\&) = Ne i{ ^-& e-^-^+^-v) , (2-2) 

where we used abbreviated notation as (7f|y) = / d 3 x Yla=o ^a( x , t)(p a (x) and 
{(p\G\ip) = J d 3 xd 3 yYlab=o ( Pa(x)G a b(x,y,t)ipf ) {y). Also, J\f is a normalization fac- 
tor. Here, the variational functions are <f a (x, t), n a (x, t), G a b(x, y, t) and S a b(x, y, t). 
In this approach, (p a represent the mean fields and G aa represent the quantum fluc- 
tuations around the mean fields <p a . The reason why we adopt the above form 
(|2-2p for the trial state with variational functions (p a (x,t), 7t a (x,t), G a b(x,y,t) and 
E ab (x,y,t) is that the (<p a (x,t),TT a {x,t)) and (G a b(x,y,t), £ ab (x,y,t)) automati- 
cally satisfy the canonical variables condition 25 ) and/or canonicity conditions 26 ) in 
terms of the time-dependent Hartree-Fock theory developed in the nuclear many- 
body theory. 

The time-dependence of the variational functions is governed by the following 
time-dependent variational principle: 

5 [ dt{V\£-H\¥) = 0, (2-3) 
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where H = f d?xTL{x). As a result, we derive a basic equation of motion for 



Sahdud^ + m 2 ab (x) - ^(pafb ) ¥b = c5 a0 



(2-4) 



where we can eliminate tt due to the honor of the canonicity condition for <p and ir. 
Here, the "mass" function m^ b (x) are defined by 



m 2 ab (x) = (m 2 + ^ 2 + ^tr S®{x,x)\ 5 ab + ^<p a {x)q> b {x) + ^(x. ./•) (2-">) 



The function S^ b \x,y) corresponds to the propagator and is formally obtained as 



S {0 \x,y) 



d 4 p 

i 



-p 2 + m 2 (x) — ie 



(2-6) 



The solution of the variational equation for G aa (x,x,t) is given in terms of the 
diagonal element of the above propagator as 



Gaa(x,x,t) = G(m 2 aa ) = S$(x,x) . 



(2-7) 



In the Hamiltonian, the direction of the spontaneous chiral symmetry breaking 
is taken in the sigma direction. Thus, the solution of the mean filed <p is assumed 
as0 

/ <Po\ 



tp = tp 



(0) 






v o y 



(2-8) 



Then, we define the sigma meson mass M and the pion mass fi, which correspond 
to the masses derived in the Hartree approximation, as 

M 2 = m 2 = ml + ^l + ^G(M 2 ) + ~G(/i 2 ) , 

M 2 = m 2 u = m 2 2 = m 2 3 = m 2 + ^ 2 + ^G(M 2 ) + ^AG(/i 2 ) . (2-9) 

From (|2-4p in the static case, we can determine the chiral condensate <po by means 
of the mean filed approximation: 



(2-10) 



Here, the first component of (p^ , (pa, in (|2-8[) is not a field operator appearing in (|2-1[) . Wc 
think that there is no confusion. 
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The solution of the above equation is given in the two phases up to the order of c: 
A [\ c 

—ifQ M 2 — y -^j^ ) (chiral broken phase) (2*11) 
tpo • (chiral symmetric phase) (2T2) 

It should be noted here that the mass [i which corresponds to the pion mass does 
not satisfy the Goldstone theorem in the chiral limit, namely, even if c = 0, then 

2.2. Pion mass and recovery of the Goldstone theorem 

To see that the Goldstone theorem is satisfied in the variational approach, we 
introduce the external source field J a coupled with the <p a : 

V! = U- S £ j J a if a . (2-13) 
a 

By this source term J a (p a , the various quantities obtained in the previous subsection 
are modified as follows: 

(p — ► <p + 5(p , 

m 2 (x) — ► m 2 (x) + 5m 2 , 

S(x,x) — ► S (0) (x,x) + SS(x,x) . (2-14) 

As a result, the equation of motion for <p is recast into the equation of motion for 

Sip as 



5ip b + Sm ab {x)p\ 



(0) 



A b + m 2 ab (x)-^xp^ 

~ (6tp a x pf ] + 4°) x 5ip b ) pf ] = J a , (2-15) 
where (pffl is a solution of Eq. (|2-4p and 



Sml b (x) = f^cpoSipo + ~tr SSJ S ab + ^ (ip a x 5<p b + Sip a x cp b ) + ^SS ab . (2-16) 

Hereafter, we consider the case that the source field J a has the following form: 

J a (x) = eSa^-^ . (2-17) 

Then, the shifted quantities Sm 2 and SS as well as Sep have the same time- and 
coordinate-dependence as J a : 

Sm 2 ab {x) = Sm 2 ab e^-^ , (2-18) 
SS ab (x,x) = Smln^tfy^-^ . (2-19) 
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Here, we have defined the polarization tensor II^\ ol \q 2 ) as 



n 



(Pol)/ 2 
ab \1 



(2-20) 



Here, let us consider the pion mass which should satisfy the Goldstone theorem 
in the chiral limit. In order to derive the pion mass in the variational approach with 
external fields, we investigate the response of the following external source field: 



For this source field, only 5<pi responds, namely, 

%/0, 8(p a = (a = 0,2,3). 
We take an ansatz for the structure of matrices 5m ab and 5S ab 

/o v o 0\ 

rj 
' 
\0 0/ 



5m 2 ab (x) 



SS ab (x) 



( (55 0\ 





\ 0/ 



(2-21) 
(2-22) 

(2-23) 



It will be understood later that this ansatz is correct. 

For a = 1, we derive the equation of motion from (|2T5p as 



where 



From (p49|) and ([S^S]) . we obtain 



5m 10 y? 



A 2 

3 1 



1 



A rr(pol) 
3^01 



(Q 2 



(2-24) 



(2-25) 



(2-26) 



Thus, from (|2-24|) and the above relation, we can derive the equation of motion for 
6(pi under the source field J%: 



3 1 



■h 



The full propagator S ab (x, y) is thus defined as 

5ip a (x)/5J b (y) = S ab (x,y) . 
Then, the mass is defined as the pole of the full propagator: 



Mi 



2,^2 

3 1 



K o1) (q) 



A 7-r(pol) 
3 JJ 01 



(0) 



(2-27) 



(2-28) 



(2-29) 
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Fig. 1. The polarization tensor is depicted diagrammatically. Here, 11^° means a part 

which is gotten rid of the external line (n) and the chiral condensate ipo (cf>o) in the above 
diagram. 



This mass corresponds to the pion mass. From Eq. (|2-20j) . the polarization tensor 
n^° l \q 2 ) appearing in Eq. (|2-27p or (|2-29[) reveals the pion and the sigma meson 
pair excitation, which is shown in FigQ] diagrammatically. 

Let us try to derive more explicit form for the pion mass M n . From (|2-7p and 
(|2-6|) . the following relation is satisfied: 

G(M 2 ) - Gifj 2 ) = [ S^I ( o \,i 5—^-5 ) (2-30) 

y ' y ' J (27r) 4 \ —p 2 + M 2 — ie -p 2 + fi 2 -iej y 1 

Thus, the polarization tensor 11^°^ can be expressed as 



r(pd)^ G{M 2 )-G(y 2 
M 2 - ix 2 

From Eqs. (|2-9p and (|2-llj) or (|2- 12|) . the following relations are derived easily: 



/i 2 = sl\^- - \(G{M 2 ) - G(fi 2 )) , (chiral broken phase) (2-32) 

V o 1V1 o 

Ac 2 A 

fj 2 = M 2 — — — - —(G(M 2 ) - G(ti 2 )) . (chiral symmetric phase) (2-33) 



Thus, we obtain 



(chiral broken phase) (2-34) 

i - ^rw) = (i - ^r } (o)) m 2 - ^ . 

(chiral symmetric phase) (2-35) 

Substituting the above relation into (12-29P and using (2-11) or (2-12) again, we finally 
obtain the pion mass in the chiral symmetry broken and chiral symmetric phases as 

fX c 

M 2 = y > (chiral broken phase) (2-36) 
Ac 2 

M 2 = M 2 j , (chiral symmetric phase) . (2-37) 
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As is expected, in the chiral broken phase, if we take the chiral limit c = 0, the pion 
mass is reduced to zero which realizes the Goldstone theorem. 18 ) As is seen in the 
polarization tensor n^°^ in Eq. (|2-20p or in FigdJ the real pion mass contains the 
effects of the pion and the sigma meson pair excitations. Also, as is seen from the 
second term of the right-hand side in Eq. (|2-29p . the polarization tensor appears in 
the denominator, so this meson pair excitation is continuously connected as 1/(1 — 



ATI, 



(pol) 
01 



/3) 



1 + Ailg ol) 



/3 + \n$° 1} /3 



\n$j°^ /3 + • • • . Thus, the meson pair 



excitation is automatically included in our approach and leads to the consistency for 
the Goldstone theorem. 

2.3. Sigma meson mass 

In this subsection, we derive the sigma meson mass as is similar to deriving the 
pion mass in the previous subsection. In order to derive the sigma meson mass, we 
need to introduce the external source field for the sigma direction 

J a (x) = edaoe^-^ , (2-38) 

and to see the response of this external field. In this case, only dipo responds: 

Stpo^O, 6<p a = (a = 1,2, 3). (2-39) 

We assume the structure of matrices as 



Sm 2 ab (x) 



fc 








°\ 





rf 














v' 





Vo 








w 



5S ab (x) 



/5S a 








^ 





5 Sir 














SSir 





V o 








6Sj 



(2-40) 



For o = 0, we obtain the equation of motion for 5(po as 



-q 2 + M 2 



Jn 



where 



5m 



oo 



X(foS(p + ^S a + ■^oS' 7r 



bm\ 2 



5ml 3 



6 ob 



From ([2 T 4T]) and ([249]), we obtain 



8m 00 tp 



l-fiT£ ol )(g2)(l + |iT n (g2)) 
where we define the reduced polarization tensor as 



• X<PqS<Po 



(2-41) 

(2-42) 
(2-43) 

(2-44) 



nn(q 2 



n 



(pol) 

11 



(Q 2 ) 



\\n^ {q 2 



(2-45) 
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(<0V/ \ {a) 



0o ) 



71 

(a) 




Fig. 2. The polarization tensor n[f and -/Ioq 01 " 1 are depicted diagrammatically in (a) and (b), 
respectively. Here, 11^° ' and -/Too '' mean parts which are gotten rid of the external line (it) 
and the chiral condensate i/pq (<^o) in the above diagrams. 



As a result, we can obtain the following relation from f)2-41|) : 

-q 2 + M + Xip 2 ■ 2 00 W '} n 6 11 W ^ 6 — ^ tyo = Jo ■ 2-46 

Thus, we can derive the sigma meson mass M a as a pole of the full propagator 
5(po/SJo: 

M 2 = M 2 + Xip 2 , ■ 2 00 ) n 6 ; 6 U ' . 2-47 

1-^(0) (l + ^n(0)) 

In this formula, the polarization tensors 77^°^ and il™ are automatically included. 

Here, 7l{^° 1 ' 1 and 77qo represent the pion-pair and the sigma meson-pair excitation, 
respectively, as is depicted in Figj2j 

§3. Pion and sigma meson masses at finite temperature 

In this section, we investigate the pion and sigma meson mass at finite temper- 
ature in this variational approach. To extend the method developed in the previous 
section to the finite temperature case, the po integral is simply replaced to the Mat- 
subara sum: 

Here, u n is the Matsubara frequency of the bosons: 

oj n = 2mrT . (n : integer) (3-2) 
Under this extension, the diagonal element of the propagator is written as 

d 3 p 



S{x,x)=iT ^ I 7^3 S ( iuJ n,p) , 
where the Fourier transform S(iuj n ,p) is given as 



(3-3) 



S(iu n ,p) = -5- 5- ^r— ~ • ( 3 ' 4 ) 

LOn + p l + m 2 (x) — ie 
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Also, the polarization tensor (12-20j) is rewritten as 

d 3 p 



n ab (iv n ,q) = -T ^2 / ( 2 ■ )3 Saa(iu m + iUn,P + q)Sbb(iu m ,p) ■ (3-5) 

m=-oo J ^ ' 



Using the well-known technique of the contour integral, the sum of the Mat- 
subara frequency is easily taken and the resultant expression for S(x, x) can simply 
expressed in terms of the bose distribution function rv^ as 



S(x,x)- — 2 ^ 



n 



(±) 



d\p\ 



1 



\Jp 2 + m 2 {x) 



n (+) _ n ( ) 



eX p ( ±V^E 



(3-6) 
(3-7) 



1 



Similarly, the polarization tensors at finite temperature are calculated and the results 
are given as follows: 



^ P a Ol) (0,0) 



p 2 + m 2 (x) 



(n (+) _ (n H) 



y/p 2 + m 2 { 



(3-8) 



n 



(pol) 
ab 



(0,0) 



Air 2 (m 2 — m 2 ) 



yjp 2 + m 2 {i 



d\p\p 2 



p 2 + m?(x) 



(-)> 



(a ^ 6) (3-9) 



where we define m a = m aa and n a = n aa for simplicity. 

§4. Numerical results and discussion 

In this section, we demonstrate the numerical results for the meson masses, M n 
and M a , and the mean filed value for (po at finite temperature. The meson masses are 
given in Eqs. (|2-29p and (|2-47p with the propagator and polarization tensors replaced 
by Eqs. (|3-3p ^ (|3-5p or (|3-6|)M|3-9|) in the finite temperature system. We have model 
parameters mo, A, c and the three momentum cutoff A. We take these value to 
reproduce the pion mass (140 MeV), the sigma meson mass (500 MeV) and the 
condensate ipo (93 MeV) at zero temperature values in each value of A. Here, the 
three momentum cutoff A is not determined in this framework. Therefore, we take 
several values of A. 

It is well known that, in the Nambu-Jona Lasinio (NJL) model, the three mo- 
mentum cutoff is taken as a rather small value around 630 MeV. 23 ) Thus, we might 
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A =650 [MeV] A =650 [MeV] 




T [MeV] T [MeV] 

(a) (b) 

Fig. 3. The sigma and pion masses are shown as functions of temperature T in (a). The order 
parameter of the chiral phase transition ipo is depicted in (b). The cutoff parameter A is taken 
as 650 MeV. 



have to adopt a rather small cutoff value. Under the above consideration, we first 
take the three momentum cutoff as 650 MeV. Thus, we neglect the polarization ten- 
sor n^°^ because this tensor represents the sigma meson pair excitation. We adopt 
the sigma meson mass at zero temperature as 500 MeV, so the pair has the energy 
1000 MeV at least. This value is beyond the cutoff parameter A. This is the reason 
why we neglect the polarization tensor i7gg in the following numerical calculation. 

Thus, we set i7gg°^ = in the formula for the sigma meson mass, Eq. (|2-47p . 

We show the pion and sigma meson masses and the mean field value ipo in Figj3] 
(a) and (b), respectively. The horizontal axis represents the temperature. In this 
case, the order parameter, ipo, of the chiral phase transition monotonically decreases 
and is a single- valued function of the temperature T, so the order of phase transition 
may be likely the second order or crossover. Correspondingly, the pion mass and 
the sigma meson mass are also monotonically changed and they are single-valued 
function of T. 

In FigU we show the results of the pion mass, sigma meson mass and the mean 
field value of the order parameter (fo for the various values of the three momentum 
cutoff parameter A. Solid curves, dash-dotted curves and the dotted curves show the 
numerical results for the cutoff parameters A = 700 MeV, 900 MeV and 1100 MeV, 
respectively. If the cutoff parameter A is taken as a rather larger value, the change of 
the order parameter is steep in the transition region. Beyond a certain value for A, 
the order parameter becomes a multi-valued function of T. Also, the meson masses 
become multi-valued functions of T. This behavior reveals the first order phase 
transition. Thus, it should be noted that the choice of the value of cutoff includes 
subtle problem such as the determination of the order of phase transition. It may 
be concluded that the rather small value for the cutoff A should be adopted in this 
model, because the chiral phase transition at finite temperature may be crossover in 
the realistic parameterization for the pion mass. 

In Figj5j the dotted curves show the results for the pion mass and sigma meson 
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A =700,900,1100 [MeV] 



A =700,900,1100 [MeV] 




Fig. 4. The sigma and pion masses are shown as functions of temperature T in (a). The order 
parameter of the chiral phase transition ipo is depicted in (b). The cutoff parameters A axe 
taken as 700 MeV (solid curves), 900 MeV (dash-dotted curves) and 1100 MeV (dotted curves), 
respectively. 



A =650 [MeV] 



A =650 [MeV] 



b 





T [MeV] 
(a) 



T [MeV] 
(b) 



Fig. 5. The sigma and pion masses are shown as functions of temperature T in (a). The order 
parameter of the chiral phase transition tpo is depicted in (b). The pion masses at zero temper- 
ature are taken as 140 MeV (solid curves) and 70 MeV (dotted curves), respectively. The cutoff 
parameter A is taken as 650 MeV. 



mass in (a) and the order parameter in (b) with the tree momentum cutoff A = 650 
MeV, where the pion mass at zero temperature is taken as 70 MeV unrealistically. 
The solid curves represent the results under the realistic value for zero temperature 
pion mass for the comparison. For lower pion mass, the order parameter becomes a 
multi- valued function of the temperature T, so the chiral phase transition becomes to 
the first order one. Thus, in the case of the chiral limit, the chiral phase transition 
is the first order one within our approach in this 0(4) linear sigma model. This 
behavior is the same as that derived in Refill), where the optimized perturbation 
theory is used in the same model. 
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§5. Summary and concluding remarks 

In this paper, we investigated the chiral phase transition at finite temperature 
in the 0(4) linear sigma model based on the time-dependent variational method 
devised by the linear response theory for the external fields. In this approach, the 
Goldstone theorem for the spontaneous chiral symmetry breaking is retained al- 
though the variational method is adopted. The reason why the above-mentioned 
situation is realized is that the meson-pair excitations are automatically included in 
this approach such as the particle-hole pair excitations for the collective motions of 
the nuclear many-body problem in the random phase approximation (RPA). It is well 
known that the pion as a Nambu-Goldstone boson has properties of the collective 
mode for the quark-antiquark pair in terms of the model based on the quarks. 27 ) 

In the 0(4) linear sigma model, the nature of the collective mode is revealed 
with the pion-sigma meson pair excitations like the RPA in our variational approach 
devised by the linear response theory. The same situation is realized for the sigma 
meson which is the chiral partner of the pion in terms of the chiral symmetry. It 
was shown that the sigma meson contains the effects of the pion-pion and sigma 
meson-sigma meson pair excitations. 

We extended this variational method to the finite temperature case by means of 
the imaginary time formalism for the finite temperature field theory. This is a natural 
extension in order to investigate the chiral phase transition at finite temperature. 

We have one free model parameter that is the three momentum cutoff parameter 
A since we regarded the 0(4) linear sigma model as a low energy effective model of 
QCD. As is similar to the NJL model, we adopted a rather small cutoff A around 700 
MeV. In this parameterization, the order parameter of the chiral phase transition 
is monotonically changed and is the single-valued function of the temperature T. 
However, the change of the order parameter in the transition region becomes rapid 
as the cutoff A is larger. Beyond a certain value of A, the order parameter becomes 
a multi-valued function of T. Thus, it should be noted that the cutoff has to be 
determined carefully. The similar caution should be given for the pion mass. It was 
shown that, even if the cutoff parameter is taken as a small value, the order of the 
chiral phase transition changes according with the pion mass at zero temperature. 
Thus, the important problem may be to determine the value of cutoff through the 
physical process, for example, low energy pion-pion scattering, if the 0(4) linear 
sigma model is used as a low energy effective model of QCD. Further, it is well 
known that a certain susceptibility presents an important information for the chiral 
phase transition. These investigations are future problems. 

Recently, the chiral phase transition at finite temperature is investigated in the 
same model by using a slight different manner, 28 ) in which the meson excitations are 
treated by introducing the Wigner functions and the Vlasov equations are formulated 
for these Wigner functions. 29 ) Thus, as another future problem, it may be important 
to understand the dynamical chiral phase transition in the several situations, for 
example, in the quench scenario, 30 - 1 in the context of the high energy heavy ion 
collisions. 
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